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We present a search for the decay t^vi {£ — r, /i, or e) in (458.9 ± 5.1) x lO*" BB pairs 

recorded with the BABAR detector at the PEP-II B-Factory. We search for these B decays in a 
sample of B^ B~ events where one B-meson is reconstructed as B~ D^l~vX. Using the method 
of Feldman and Cousins, we obtain S(B+ ^ r+v-r) = (1.7 ± 0.8 ± 0.2) x lO""*, which excludes zero at 
2.3a. We interpret the central value in the context of the Standard Model and find the B meson decay 
constant to be /| = (62 ± 31) x 10^ MeV^ We find no evidence for 5+ e+i/e and B+ fi'^iy^ 
and set upper limits at the 90% C.L. B{B+ -> e+i/g) < 0.8 x lO"'^ and B{B+ ^ fJ-^i^^) < 1.1 x 10"^ 

PACS numbers: 13.20.-v, 13.25.Hw, 12.15.Ji 

In the Standard Model (SM), the purely leptonic decay £~^i>e f]\ proceeds via quark annihilation into a 
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boson. This process is related to the Cabibbo-Kobayashi- 
Maskawa matrix element Vub and the B meson decay 
constant, Jb, by B{B+ — > i^i^e) |Ki6p/s- It is also 
potentially sensitive to the presence of a charged Higgs 
boson as in the minimal supersymmetric extension 
of the Standard Model. Using \Vub\ = (3.94 ± 0.26) x 
10-3 0] and fs = 190 ± 13 MeV [J and assuming 
only a SM contribution to the process, the branching 
fraction predictions are B{B^ t^i^t) = (1-0 ± 0.2) x 
10-^ B{B+ fi+j^^) = (4.5 ±1.0) X 10-^ and B{B+ 



-Ve) = (1-1 ± 0.2) X 10 
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The different branching 



fractions result from helicity suppression of the lower- 
mass charged leptons. In this paper, we describe a search 
for all three final states. 

The data used in this analysis were collected with the 
BaB^ detector at the PEP-II storage ring at the SLAC 
National Accelerator Laboratory. We use the full BaBAR 
dataset, corresponding to an integrated luminosity of 
417.6 fb""^ with center-of-mass (CM) energy equal to the 
r(45) rest mass. These data contain (458.9 ±5.1) x 10^ 
T{AS) BB pairs and we assume equal production of 
BW and B+B- from the T{4S) decays. The BaBAR de- 
tector is described in detail elsewhere [Hj. For the most 
recent 203 fb^^ of data, the barrel region of the muon 
system was upgraded to limited streamer tubes Q . 

Signal and background processes are simulated us- 
ing EVTGEN 0. A GEANT4-based Monte Carlo 
(MC) simulation is used to model the detector response 
and to estimate the signal efficiency and the physics 
backgrounds. Simulation samples equivalent to approxi- 
mately three times the accumulated data were used to 
model BB events, and samples equivalent to approx- 
imately 1.5 times the accumulated data were used to 
model continuum background events where e'^e~ — > uu, 
dd, ss, cc, and t~^t~. We independently simulate the 
signal processes at a rate over a hundred times that ex- 
pected in data, using samples where one B meson always 
decays as B^ t^vi and the second decays into any 
final state. We normalize these signal samples to their 
predicted SM branching fractions. 

The strategy adopted for this analysis is similar to that 
from our previously published work [9|. Signal B decays, 
_B+ i^i'e, are selected in the recoil of a semileptonic 
decay, B' D^i-vX, referred to as the "tag" B. The 
final states of the r+ decay in i?+ r^Vr are identical 
to those in Ref. [9||: r+ ^ e'^VeT'T, r+ /.t+f^iF,-, r+ — > 
tt^^Ft, and t+ — > p^Vr- For the first time, we include 
5+ — * e+j/g and 5+ i-^-^v^ in this search. In addition 
to using about 20% more data than in Ref. [9| , we relax 
the constraints on the tag B, improve the definition of 
the discriminating variables and use a combination of tag 
and signal B variables in a multivariate discriminant that 
improves signal efficiency and background rejection. 

The tag B is reconstructed in the set of semileptonic B 
decay modes B~ ^ D^i~i'X, through the full hadronic 
reconstruction of D'^ mesons and identification of the lep- 



ton, £~ , as either e~ or ^~ . Other particles {X) resulting 
from a transition from a higher-mass charm state down 
to the D'^ are not explicitly reconstructed and are not 
included in the tag B kinematics. This strategy, and 
the reconstruction method (Z)" decay modes, D'^l^ ver- 
tex requirements, etc.), are the same as in Ref. Q. One 
difference in the present analysis is that we may assign 
up to one photon (from X) back to the tag B, based on 
its consistency with the decay D*° — > (7r°, 7)!?°. 

The efficiency for tag B reconstruction (etag) is defined 
as the rate at which events in the signal MC are found 
to contain at least one reconstructed tag B and a single 
track recoiling against that tag. The efficiency for each 
signal mode is given in Table llllj including corrections 
for systematic effects (described below). The efficiency 
is larger for _B+ —^ t^Vt events due to high- multiplicity 
T+ decays faking tag B mesons. 

We identify one of the following reconstructed parti- 
cles recoiling against the tag B: e+, 7r+, or yO+. The 
e"*" and /x"*" can come from B~^ t^i^t, with the r"*" 
decaying leptonically, or directly from B^ l^-^^tj, or 
B"*" e^Ve- The signal track must originate from the 
interaction point (IP), with a distance of closest approach 
to the IP less than 2.5 cm along the beam axis and less 
than 1.5 cm transverse to the beam axis. We reject 
events that contain more than one such IP track recoil- 
ing against the tag B. There may be additional tracks 
that do not come from the IP. We reject events where 
the single IP track is identified as a kaon. We assign 
the single-track recoils to categories based on a hierar- 
chical selection. An event is assigned to the /i+ cate- 
gory if the track passes muon identification or to the e+ 
category if it passes electron identification; in the latter 
category, we recover up to one bremsstrahlung photon 
based on angular separation from the track and add its 
four-momentum to the electron's. We assign the event 
to the category if it fails lepton identification and 
can be paired with a tt*^ candidate. The tt*^ candidates 
used in the p"*" reconstruction are defined as a pair of 
photons, each with laboratory energy > 50 MeV, with 
invariant mass m^a — [0.115, 0.150] GeV/c^. Single-track 
events that fail the selections above are assigned to the 
TT"*" category. 

While the direction of neither B meson can be known 
precisely, four-momentum conservation constrains the 
tag B momentum to lie on a cone around the flight di- 
rection of the reconstructed D^i" system. The cosine of 
the opening angle between the B meson and the D'^i^ 
system in the CM frame is given by 



cos Ub y 



IEbEy 



2|PB||pi 



(1) 



where Y refers to the reconstructed tag B final state, 
[Ey, Py) and {Eb, Pb) are the four-momenta in the CM 
frame, and my and uib are the masses of the Y system 
and tag meson, respectively. Eb and the magnitude 
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oipB are calculated from the beam energy: Eb = -Ecm/2 
and \pb\ = \/ Eg — rrig . Decays of the B meson directly 
to D^i~v are largely constrained to the physical region of 
this cosine, while decays involving a higher-mass charm 
state will yield cosine values below the physical region 
when the intermediate decay particles (e.g. tt*^ or 7) are 
not explicitly reconstructed. 

The signal B momentum vector is equal in magnitude 
to \pb\ and is opposite to the tag B direction, so that 
it lies on the cone of the tag B momentum defined by 
Equation [TJ To estimate quantities in the signal B rest 
frame, such as the momentum of the signal B daugh- 
ter (s), we choose the signal B boost vector on that cone 
and we compute the quantity in the corresponding rest 
frame. We then use the value of that quantity averaged 
over all trial rest frames as an estimate of the true value. 
We denote the momentum of the signal particle(s) de- 
termined by this method as p^jg. This has the largest 
impact in the 5+ e^i^e and _B+ l^^v^ channels, 
where the Icpton is mono-energetic in the signal B rest 
frame. The improved resolution of the lepton momen- 
tum directly improves the separation of signal and back- 
ground. If an event has a reconstructed signal muon 
(electron) candidate and Pgjg > 2.30 (2.25) GeV/c, it is 
classified as a 5+ f^^i^t^ [B^ e+i^e) candidate; oth- 
erwise, it is classified as _B+ t^i^t, with t+ fi'^VfjVr 

(t+ e'^VeVr). 

A critical discriminating variable is the extra energy 
(i?oxtra), which is the total energy of charged and neutral 
particles that cannot be directly associated with the re- 
constructed daughters of the tag B or the signal B. This 
variable was not examined (kept "blind" ) until the analy- 
sis strategy was finalized. We expect the signal to concen- 
trate near zero -EcxtraJ however, due to collider-induced 
backgrounds, detector noise, and unassigned tracks and 
neutrals from the tag and signal B mesons, signal events 
can have non-zero -Ecxtra- We require a minimum energy 
in the laboratory frame of 30 MeV for any neutral cluster 
used in i?extra- We improve our signal and background 
separation in this variable by using an algorithm to as- 
sign up to one photon from the E'cxtia back to the tag B. 
Candidate extra photons must have a CM-frame energy 
less than 300 MeV, consistent with having come from a 
tt" or 7 from the D*^ transition. If, by adding a 

candidate photon back to the tag B kinematics, the value 
of cos y becomes closer to (but not greater than) 1.0, 
it is retained as a transition particle candidate. If more 
than one photon satisfies these conditions, the one which 
moves AM = rnno-y — mjjo closest to the nominal value 
of 142 MeV/c^ [10| is used. This photon is excluded from 
Ecxtra- The tag B kinematic quantities and p'^-^^ are re- 
computed, with the photon added to the tag B final state. 

The background consists primarily of B^ B~ events in 
which the tag B meson has been correctly reconstructed 
and the recoil contains one reconstructed track and ad- 



ditional particles that are not reconstructed. Typically 
these events contain mesons and other particles that 
are not detected and thus fake the multiple neutrinos in 
signal events. Backgrounds from B decays and contin- 
uum processes have distinctive signatures in a number of 
discriminating quantities. We group variables according 
to those which are computed from the whole event, the 
tag B, the signal B, and other sources. Some variables, 
such as those associated with the whole event, are useful 
for rejecting continuum background, while others (such 
as those associated with the reconstructed B mesons) are 
better at rejecting B background. 

The event-level variables are: the ratio of the second 



and zeroth Fox- Wolfram moments ll| ; the minimum in- 
variant mass of any two charged tracks in the event; the 
net charge of the event; cos9b,y', the invariant mass of 
the two leptons in the event (mu); and the missing mass 
vs. cosine of the polar angle (laboratory frame) of the 
missing three-momentum, where the sum defining the 
reconstructed four-momentum runs over all charged and 
neutral particles in the event. The tag B variables are: 
the D'^ decay mode; the CM momenta of the tag B kaon 
and lepton; particle identification quality of the tag B 
charged kaon (where applicable). The signal B variables 
are: the quality of the particle identification of the signal 
muon, for muon final states of the signal B; the quality of 
the kaon identification on the signal track (to reject kaons 
misidentified as leptons or pions); for r"*" — s- 7r"*'7r'^I7T-, 
the reconstructed mass of the , and the CM momenta 
of the /9+ daughters; and for _B+ t^i't, cosO'^y 
Pgj^, where cos 6*^ y defined in the signal B meson rest 
frame using Equation [U replacing B meson quantities 
with those of the t {Er = mB/2 and Pr = — ml) 

and where Y refers to the reconstructed r final state 
(computed using the signal B meson rest frame averag- 
ing procedure). Other variables used are: the separation 
between the tag B meson decay vertex and the point of 
closest approach to the IP by the signal B track; and the 
distribution of the cosine of the angle between the signal 
B CM momentum and the tag B thrust vs. the mini- 
mum invariant mass of any three charged particles in the 
event Q. 

The shapes of these variables in MC simulation are 
then used to define probability density functions (PDFs) 
for signal (Pg) and background [Pb)- We define for each 
variable the ratio Ps/[Pb + Ps\- We use the product 
of these ratios to construct a pair of likelihood ratios 
(LHRs) for each signal channel, one for rejecting B back- 
grounds {LHRg-g) and the other for rejecting continuum 
{LH Rcont.) backgrounds. The LHR output is bounded 
between and 1, with signal accumulating toward 1 and 
background toward 0. 

We optimize selection criteria on Scxtra, LHRg-g, and 

> e^Vf, and 



LHRc 

B+ - 



for all modes. For the B^ 



modes, we additionally optimize the se- 



lection on Pgjg. 



For the r+ 



mode we addi- 
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tionally optimize the selection on m« (to reject poorly 
modeled photon-conversion background). For the t de- 
cay modes, we choose the figure-of-merit (FOM) to be 
^sig/ \/ {Nsig + -/Vbg), since there is still significant back- 
ground left in these channels even after final selection 
criteria are applied. For i?+ — > /i^i'^ and — > e+z^e 
we use A'^gig/(3/2 -|- ^/N\^) [1^] due to the low expected 
background. We divide the MC simulation samples for 
signal and background into thirds, two for optimization 
and one from which to compute unbiased efRciencies and 
background predictions. This latter sample has statis- 
tics roughly equivalent to the data. Optimized selection 
criteria are given in Table [H The signal efficiency (ffgig) 
is defined as the rate at which signal events containing 
a reconstructed tag B are also found to contain a signal 
B candidate, and it includes the r"*" branching fractions. 
These efficiencies are given in Table IIIIl 



TABLE I: Optimized signal selection criteria. 



Mode 


LHRgg 


LH Rcont. 














(GeV) 


( GeV/c) 


(GeV/c2) 






B+ - 


-> T~^U-r 






e'^uu 


> 0.77 


> 0.25 


< 0.20 




> 0.29 




> 0.14 


> 0.72 


< 0.24 








> 0.97 


> 0.95 


< 0.24 








> 0.57 


> 0.80 


< 0.35 










B+ ^ 


{tJ.+ ,e+)u 






> 0.33 


> 0.61 


< 0.72 


[2.45,2.98] 






None 


> 0.01 


< 0.57 


[2.52,3.02] 





We calibrate our background prediction using side- 
band regions of Scxtra where the signal contribution is 
negligible. We define the sidebands for 5+ T'^Vr, 
B+ fJ.'^i^fj., and B+ e+z/g as E'extra > 0.4 GeV, 
> 0.72 GeV, and > 0.6 GeV, respectively. We predict 
TV'i^g*'', the number of background events in data in the 
E'oxtra signal region (Table HI)), by scaling the yield pre- 
dicted by the MC simulation {N^'^) by the ratio of yields 
in data {N^X,"") and MC (iV^g) in the sideband. This 
method assumes that the shape of E'extra is well-described 
but does not rely on the absolute prediction of the yield. 
We validate this approach by defining sidebands in other 
variables (D" mass, LH Rcont., LHRg-g, and Psjg) and 
studying the data/MC agreement for the entire Soxtra 
background shape. We find the shape to be well de- 
scribed. We also studied the effect of varying the i?oxtra 
sideband definition and obtained consistent background 
predictions. 

The branching fraction for any of the decay modes is 



obs 



A /"data 
■■'^bg 



2Nb + B-£tag£sig 



(2) 



where iVobs is the total number of events observed in the 
signal region and Nb+b- is the total number of T{AS) 



TABLE II: Background predictions from the iSextra sideband, 
as described in the text. 



Mode 


^ ^sidc 


T\ 7data 




* /"data 


r+ -> 


e'^VeVr 


333 ± 19 


334 ± 18 


81 ± 10 


81 ± 12 


T+ -> 




1248 ± 36 


1236 ± 35 


136 ± 12 


135 ± 13 


r+ ^ 




6507 ± 88 


7167 ± 85 


212 ± 19 


234 ± 19 


r+ ^ 


P+Vr 


1841 ± 48 


1734 ± 42 


62 ± 9 


59 ± 9 


B+~. 




12 ± 5 


14 ± 4 


12 ± 5 


13 ± 8 


B+ ^ 




26 ± 6 


42 ± 6 


15 ± 5 


24 ± 11 



B^B decays in the data. The estimation of Nb+b- has 
an uncertainty of 1.1% [l3l |. 

Potential sources of significant systematic uncertainty 
in £tag and egig include the tag reconstruction rate, the 
modeling of i?oxtra, and signal track and neutral recon- 
struction. We use "double-tagged" events to study possi- 
ble effects. Double-tagged events contain two fully recon- 
structed, independent, oppositely charged semileptonic 
tag B decays. These double-tagged events are analogous 
to signal, in that every particle that can be assigned to 
the original B decays has been assigned. 

We use the absolute yields of tagged events to obtain 
a systematic uncertainty on etag- We form a double ratio 
from the ratios of double-tagged to single-tagged events 
in the data and MC simulation. Single-tagged events 
are defined as events containing at least one semilep- 
tonic tag B decay with no constraints on the rest of the 
event. We improve the sample purity by requiring that 
— > K^n^ in at least one of the tags. We measure this 
double ratio to be 0.891 ± 0.021. As a comparison, we 
perform the same measurement replacing K~t:'^ 
with /C~7r+7r~7r+ and find the double-ratio to be 

0.954 ± 0.011. We use 0.891 as the nominal correction 
to Stag and treat the relative difference between the two 
methods (7.1%) as the systematic uncertainty. 

The i^cxtra distribution in double-tag events is ex- 
pected to contain contributions similar, though not iden- 
tical, to those from signal events. We validate £^extra us- 
ing the double-tagged events described above, addition- 
ally requiring that the second tag contains only — > 
K^TT^ and satisfies cos6b,y = [—1-1, 1-1] to reject second 
tags with missing neutrals. The resulting i?cxtra distri- 
bution is shown in Fig. [T] It is well-described by the MC 
simulation. We compare the efficiency of selecting events 
in data and MC simulation for i?cxtra < 0.4 GeV and find 
that the efficiency needs to be corrected by 0.985 ± 0.044 
to match the data. The uncertainty on this correction 
is due to the statistical uncertainty on the data and MC 
simulation, and we treat it as a systematic uncertainty. 

The remaining systematic uncertainties on Sgig come 
from tracking efficiency (0.36% per signal track), tt" re- 
construction for the r+ — > p'^T'r mode (0.984 ± 0.030), 
and particle identification. These are evaluated using 
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Ej„„ from Double Tag Sample (GeV) 

FIG. 1: Distribution of iScxtra in double-tagged events. The 
data (black points) and MC simulated events (gray rectan- 
gles) are normalized to unit area. The rectangles represent 
the MC simulation uncertainty. 



control samples of well-characterized particles. The par- 
ticle identification efficiency corrections and systematic 
uncertainties are 0.953 ±0.003 (0.97 ±0.04) for identified 
electrons in the 5+ t^Vt {B~^ e+z/g) analysis and 
0.92 ± 0.05 (1.016 ± 0.022) for identified muons in the 
5+ t^Vt (5+ l^^i^fj.) analysis. 
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TABLE III: The corrected tag and signal efficiencies. The first 
uncertainty is the MC statistical uncertainty, and the second 
is the systematic uncertainty from sources described in the 
text. Branching fractions are included (e.g. — > e'^uV). 
The last column is the total systematic uncertainty on each 
efficiency as a percent of its value. 



Channel 


Efficiency (%) 


Uncertainty (%) 


Tag Efficiencies 


B+ T+Vr 


(1.514 ±0.003 ±0.107) 


7.1 


B+ fi+u^ 


(0.937 ± 0.003 ± 0.066) 


7.1 


B+ ^ e+z/e 


(0.974 ± 0.003 ± 0.069) 


7.1 


Signal Efficiencies 


— + e^Vel'r 


(1.58 ±0.04 ±0.07) 


4.5 




(1.45 ±0.03 ±0.11) 


7.4 


r"*" —^ -K^Ut 


(2.44 ±0.05 ±0.11) 


4.5 


r+ p+Vr 


(0.83 ±0.03 ±0.05) 


5.4 


B+ T+Vr 


(6.31 ±0.07 ±0.34) 


5.4 


B+ fi+u^ 


(28.65 ± 0.34 ± 1.75) 


6.1 


B+ e+u. 


(37.01 ± 0.38 ± 1.84) 


5.0 



The Soxtra distributions for each channel are given in 
Fig. [2] and results given in Table HVl We use the method 
of Feldman and Cousins [3] to interpret the yields in 
each channel. When computing the level at which we 
exclude the null hypothesis, we include systematic er- 
rors as a Gaussian convolution with the nominal Pois- 
son distribution. Our results in the B+ fJ^+v^ and 



FIG. 2; iScxtra after all selection criteria have been applied 
for each final state. Shown are data (black points), back- 
ground MC simulation (gray shaded), and signal MC sim- 
ulation (dotted line) normalized to 10 times the expected 
branching fraction (10^ times for B+ — > e+Ue). The back- 
ground MC simulation is luminosity normalized and corrected 
for the data/MC ratio in the iScxtra sideband; the rectangles 
represent the MC simulation statistical uncertainty. In (a-d), 
the vertical dashed line indicates the signal region boundary. 
In (f-g) the first bin is the signal region. 



B+ — > e+Ve channels are consistent with the background 
expectation and we obtain only one-sided 90% confidence 
intervals. For B+ t+Vt, we obtain a two-sided 68% 
confidence interval and exclude the null hypothesis at 
the level of 2.3(7. This result supersedes that of the pre- 
vious work 0. The statistical consistency test of the 
results over the four B+ — > t+Vt channels has a per 
degree-of-freedom of 2.02/3, or a probability of 57%, and 
is performed using branching fractions computed with 
Equation [21 In the context of the SM we determine that 
fg = (62 ± 31) X 10"^ MeV^, where the uncertainty arises 
dominantly from this measurement and | Vub \ ■ 

We obtain a single BABMi result for B+ r+Vr by 
combining this result with B{B+ t+Vt) = (l.Slg g) x 
10~^, which is derived from a statistically-independent 
sample usin g ta g B mesons decaying into fully hadronic 



final states [15|. We use a simple error- weighted aver- 
age, since the correlated systematics (mainly due to par- 
ticle identification, charged particle tracking, and i?cxtra) 
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TABLE IV: The expected background, observed events in 
data, and branching fraction results, determined as described 
in the text. 



Mode 


J^bg 


A^obs 


Branching 






fraction (xlO~*) 


r"*" — > e^VeVr 


81 ± 12 


121 


(3.6 ± 1.4) 




135 ± 13 


148 


r+ p+Fr 


59 ± 9 


71 


(2.1^?:°) 


r+ ^ ^+-r 


2.34 ± 19 


243 




B ~ t' jy^ 


509 ± 30 


583 


(i.7 ^ 0.8 ± 0.2) 


B+ ^ 11+ 1^^ 


13 ± 8 


12 


< 0.11 (90% CL.) 


B+ e+Ve 


24 ± 11 


17 


< 0.08 (90% CL.) 



have a negligible impact on the combination. Wc obtain 
B{B+ T+Ur) = (1.7 ± 0.6) X 10-*, which excludes zero 
at the 2.8cr level. Both this and the combined results are 
consistent with the SM prediction. 

In conclusion, we have used the complete BaBAR data 
sample to search for the purely leptonic B meson decay 
5+ l^v using a scmileptonic B decay tagging tech- 
nique. We measure B{B+ r+Vr) = (1.7 ± 0.8 ± 0.2) x 
10~^ and exclude the null hypothesis at the level of 2.3cr. 
We find results consistent with the background predic- 
tions for the decays — > iX^v^ and B^ — > e+i^e. We are 
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tions that support BaBAR. The collaborating institutions 
wish to thank SLAC for its support and kind hospitality. 
This work is supported by DOE and NSF (USA), NSERC 
(Canada), CEA and CNRS-IN2P3 (Prance), BMBF and 
DFG (Germany), INFN (Italy), FOM (The Netherlands), 
NFR (Norway), MES (Russia), MEG (Spain), and STFG 
(United Kingdom). Individuals have received support 
from the Marie Guric EIF (European Union) and the 
A. P. Sloan Foundation. 



' Now at Temple University, Philadelphia, Pennsylvania 

19122, USA 

* Also with Universita di Perugia, Dipartimento di Fisica, 
Perugia, Italy 

^ Also with Universita di Roma La Sapienza, 1-00185 
Roma, Italy 

' Now at University of South Alabama, Mobile, Alabama 
36688, USA 

** Also with Laboratoirc dc Physique Nuclcairo ct do 
Hautes Energies, IN2P3/CNRS, Univorsitc Pierre et 
Marie Curie-Paris6, Universite Denis Diderot-Paris7, F- 
75252 Paris, Prance 

Also with Universita di Sassari, Sassari, Italy 
[1] Charge-conjugate modes are implied throughout this pa- 
per. 

[2] W.-S. Hou, Phys Rev. D48, 2342 (1993). 

[3] E. Barberio et al. (Heavy Flavor Averaging Group), 

arXiv:0808.1297 [hep-ex] (2008). 
[4] E. Gamiz et al. (HPQCD Collaboration), Phys. Rev. 

D80, 014503 (2009). 
[5] B. Aubert et al. (BABAR Collaboration), Nucl. Instrum. 

Methods A479, 1 (2002). 
[6] W. Menges, IEEE Nucl. Sci. Symp. Conf. Rec. 5, 1470 

(2006). 

[7] D. J. Lange, Nucl. Instrum. Methods A462, 152 (2001). 

[8] S. AgostinelU et al. (GEANT4 Collaboration), Nucl. In- 
strum. Methods A506, 250 (2003). 

[9] B. Aubert et al. {BABm Collaboration), Phys. Rev. D76, 
052002 (2007). 

[10] C. Amsler et al. (Particle Data Group), Phys. Lett. 

B667, 1 (2008). 
[11] G. C. Fox and S. Wolfram, Phys. Rev. Lett. 41, 1581 

(1978). 

[12] G. Punzi, Proceedings of PHYSTAT2003: Statistical 
Problems in Particle Physics, Astrophysics, and Cosmol- 
ogy pp. 79-83 (2003). 

[13] B. Aubert et al. {BABAR Collaboration), Phys. Rev. D67, 
032002 (2003). 

[14] G. J. Feldman and R. D. Cousins, Phys. Rev. D57, 3873 
(1998). 

[15] B. Aubert et al. {BABm, Collaboration), Phys. Rev. D77, 
011107 (2008). 



Deceased 



